home *** CD-ROM | disk | FTP | other *** search
/ NetNews Offline 2 / NetNews Offline Volume 2.iso / news / comp / lang / c-part1 / 8570 < prev    next >
Encoding:
Text File  |  1996-08-05  |  3.0 KB  |  115 lines
  1. Path: keats.ugrad.cs.ubc.ca!not-for-mail
  2. From: c2a192@ugrad.cs.ubc.ca (Kazimir Kylheku)
  3. Newsgroups: comp.lang.c
  4. Subject: Re: Fortran to C conversion help?
  5. Date: 4 Mar 1996 19:15:03 -0800
  6. Organization: Computer Science, University of B.C., Vancouver, B.C., Canada
  7. Message-ID: <4hgbjnINN136@keats.ugrad.cs.ubc.ca>
  8. References: <4hf04k$ogs@pipe3.nyc.pipeline.com>
  9. NNTP-Posting-Host: keats.ugrad.cs.ubc.ca
  10.  
  11. In article <4hf04k$ogs@pipe3.nyc.pipeline.com>,
  12. Doug Kolmar <dkolmar@grovestocktn.com> wrote:
  13. > 
  14. >Hello folks, 
  15. > 
  16. >I'm sort of a weekend C programmer so perhaps I've missed something, but
  17. >I'm trying to convert a function in FORTRAN to C. 
  18. > 
  19. >Attached is the FORTRAN code followed by the C code that I've written. The
  20. >C program will compile and run, but does not yield the correct output. I.E.
  21. >integers in the range 0 to 31. 
  22. > 
  23. >Except for the differences in the range of the RAND functions, I know
  24. >nothing about FORTRAN, so am I missing something implicit in the code? Any
  25. >help would be greatly appreciated. 
  26.  
  27. Use an automatic translator. Here is the output of f2c:
  28.  
  29. First of all, I had to fix up your fortran program. Fortran requires that the
  30. first six characters of each punched card^H^H^H^H^H^H^H^H^H^H^H^Hsource line
  31. be blank, or contain a line number or 'C' character for comment. Secondly, you
  32. are missing a line 20, the target of the GO TO statement. I just stuck in a
  33. line 20 wherever I liked:
  34.  
  35. file fortran.f:
  36.  
  37.  
  38.       FUNCTION I1OVRF(LAST) 
  39. C     THE ARGUMENT, LAST, IS THE PREVIOUS VALUE OF THE SEQUENCE 
  40.       NEW=0 
  41. C     THE VARIABLE K EQUALS ONE-HALF THE NUMBER OF POSSIBLE VALUES 
  42.       K=16 
  43.       L=LAST 
  44. C     THE VARIABLE PROBIT EQUALS 1/(NUMBER OF POSSIBLE VALUES) 
  45.       PROBIT = .03125 
  46. 20    J=L/K 
  47.       IF (J.EQ.1) L=L-K 
  48.       U=RAND(0) 
  49.       IF (U.LT.PROBIT) J=1-J 
  50.       NEW = NEW+J*K 
  51.       K=K/2 
  52.       PROBIT=PROBIT*2 
  53.       IF (K.GE.1) GO TO 20 
  54.       I1OVRF=NEW 
  55.       RETURN 
  56.       END  
  57.  
  58. file fortran.c produced by AT&T's f2c:
  59.  
  60.  
  61. /* fortran.f -- translated by f2c (version 19940305).
  62.    You must link the resulting object file with the libraries:
  63.     -lf2c -lm   (in that order)
  64. */
  65.  
  66. #include "f2c.h"
  67.  
  68. /* Table of constant values */
  69.  
  70. static integer c__0 = 0;
  71.  
  72. integer i1ovrf_(last)
  73. integer *last;
  74. {
  75.     /* System generated locals */
  76.     integer ret_val;
  77.  
  78.     /* Local variables */
  79.     extern doublereal rand_();
  80.     static integer j, k, l;
  81.     static real u, probit;
  82.     static integer new_;
  83.  
  84. /*     THE ARGUMENT, LAST, IS THE PREVIOUS VALUE OF THE SEQUENCE */
  85.     new_ = 0;
  86. /*     THE VARIABLE K EQUALS ONE-HALF THE NUMBER OF POSSIBLE VALUES */
  87.     k = 16;
  88.     l = *last;
  89. /*     THE VARIABLE PROBIT EQUALS 1/(NUMBER OF POSSIBLE VALUES) */
  90.     probit = (float).03125;
  91. L20:
  92.     j = l / k;
  93.     if (j == 1) {
  94.     l -= k;
  95.     }
  96.     u = rand_(&c__0);
  97.     if (u < probit) {
  98.     j = 1 - j;
  99.     }
  100.     new_ += j * k;
  101.     k /= 2;
  102.     probit *= 2;
  103.     if (k >= 1) {
  104.     goto L20;
  105.     }
  106.     ret_val = new_;
  107.     return ret_val;
  108. } /* i1ovrf_ */
  109.  
  110.  
  111. This is one of those rare occurances when the output of a translator is more
  112. readable than the source.
  113. -- 
  114.  
  115.